Controllable Coupling in Phase-Coupled Flux Qubits 
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We propose a scheme for tunable coupling of phase-coupled flux qubits. The phase-coupling 
scheme can provide a strong coupling strength of the order of Josephson coupling energy of Josephson 
junctions in the connecting loop, while the previously studied inductive coupling scheme cannot 
provide due to small mutual inductance and induced currents. We show that, in order to control 
the coupling, we need two dc-SQUID's in the connecting loop and the control fluxes threading the 
dc-SQUID's must be in opposite directions. The coupling strength is analytically calculated as a 
function of the control flux at the co-resonance point. 

PACS numbers: 74.50.+r, 85.25.Am, 85.25. Cp 



I. INTRODUCTION 



Superconducting Josephson junction qubit is one of 
the most promising candidates for implementing quan- 
tum computation!. Single qubit coherent oscilla- 
tions in superconducting qubits have been demonstrated 
experimentally 2,3 ^ and furthermore two qubit coupling 
and entanglement have been performed in charge^, 
g U3i 7 i 8 i 9 anc j p nasQ io qubits. Scalable quantum com- 
puting requires controllable and selective coupling be- 
tween two remote as well as nearest neighbor qubits. Re- 
cently much theoretical efforts have been devoted on the 
study about the controllable coupling of charge 1 !*, charge- 
phased and flux qubits^i^i^i^. For flux qubits the con- 
trollable coupling schemes use inductive coupling, but it 
is too weak to perform efficient two-qubit gate operations. 
Hence, while in superconducting charge qubit two-qubit 
coherent oscillations and CNOT gate operations were 
experimentally observed^, only spectroscopy measure- 
ment was done for inductively coupled flux qubiA In 
this study thus we suggest a scheme to give both strong 
and tunable coupling between two phase-coupled flux 
qubits. The phase-coupling scheme, which we previously 
proposed^, has been realized in a recent experiment* 1 ^ 
The controllable coupling scheme using phase-coupled 
qubits with threading AC magnetic field was also studied 
theoreticallyiS. Further, there have been studies about 
somewhat different phase-coupling schemes22i2i. 

Two current states of a flux qubit are characterized 
by the induced loop current related with the phase dif- 
ferences across Josephson junctions in the qubit loop. 
If we try to couple two flux qubits using mutual induc- 
tance, the coupling strength J = MIlIr k 0.5GHaS will 
be too weak to perform the discriminating CNOT gate 
operations^, since the mutual inductance M and the in- 
duced currents of the left (right) qubit Il(R) 1S verv small. 
Even though the induced currents of flux qubits are weak, 
the phase differences across Josephson junctions <f> are 
as large as <f>/2-K > 0.16. Hence, if two flux qubits are 
coupled by the phase differences between two Josephson 
junctions of different qubits, we can achieve a strong cou- 
pling of the order of Josephson coupling energy E'j of the 
Josephson junctions in the connecting loop whose typical 



value is as large as up to about 200GHz. 

Introducing two dc-SQUID's interrupting the connect- 
ing loop as shown in Fig. ^we can control the coupling 
between phase-coupled flux qubits. The control fluxes, 
f' L and f'n, threading two dc-SQUID loops must be in 
opposite directions in order to give rise to the control- 
lable coupling. When two fluxes are in the same direc- 
tion, the change of control fluxes induces an additional 
current flowing in the connecting loop, causing the shift 
of qubit states as well as the change of coupling strength. 
Such a dilemma also persists in the case of one dc-SQUID 
loop in connecting loop. However, if the control fluxes are 
in opposite directions, we have found that the additional 
currents coming from two dc-SQUID's are cancelled each 
other and thus the coupling strength can be tunable re- 
maining the qubit states unchanged. 



II. PHASE-COUPLING OF FLUX QUBITS 

The three- Josephson junctions qubits 22 23-24 in Fig. 
^ has two current states; if the qubit current I cx 
—Eji sin <pi < 0, it is diamagnetic while, if / > 0, param- 
agnetic. Introducing the notation | |) (| ])) for diamag- 
netic (paramagnetic) current state of a qubit in pseudo 
spin language, there can be four current states of coupled 
qubits, | ||), | |T), I IT) and | U), of which we show one 
of the same current states, | ||), and one of the different 
current states, | ||), in Fig. ^ The phase 4>li and 4>m of 
the Josephson junctions of the three- Josephson junctions 
qubits have different values if two qubits are in different 
states. Then the phase difference 4>li — 4>Ri induces the 
phases ^ in the Josephson junctions of dc-SQUID loops. 

If we neglect small kinetic inductance, the boundary 
conditions of the left (right) qubit and the connecting 
loop can be approximately written as 

<t>L(R)l + <f>L(R)2 + 0L(B)3 = 27r(n i(fl) + f t ,L(R)), (1) 

# + & = 2n(r + f[ nd ) + (0 L1 - fa), (2) 
-J>[ + (t>' 2 = 2ir(f L + P ), -& + & = -27r(/fl + 9), (3) 



where f tL 



(R) 



L(R) 



find 



L(R) 



is total flux and 



h(R) = $ext,£(n)/$o with the external flux $ e xt,z(ii) 
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and the unit flux quantum $o = h/2e is dimension- 
less reduced flux threading the left (right) qubit. Here 
/ind,L(ii) = L s Il(r)/$o with the self inductance L s and 
the induced current Il(r) of qubit loop is the induced flux 
of each qubit and // nd = L' S V /$o that of the connecting 
loop and nL,nn, r,p and q are integers. 

We consider that the external fluxes /x and /jj thread- 
ing the qubit loops are also in opposite directions, since 
they are connected in a twisted way in the scalable de- 
sign of Ref. 0. However, for just two qubit coupling, 
we can choose the directions of external fluxes threading 
the qubit loops arbitrary. Actually there is no exter- 
nal flux in the connecting loop, but the phase difference 
{4>li — 4>ri) in the boundary condition of Eq. J5J plays 
the role of effective flux in the connecting loop, 



fefi 



?L1 



*>R1 



2tt 



(4) 



When two qubits are in different current state, i.e., one is 
diamagnetic and the other paramagnetic, the value of f' eS 
becomes 0.3 ~ 0.7. Since the induced flux of flux qubit 
is so weak as / ind rj 0.002, large value of f' eS » /j nc i in 
the phase-coupled flux qubits can give a strong coupling 
compared to the inductive coupling scheme. 

The Hamiltonian of the coupled qubits can be given 

by 
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H = -P J ■ M • P + E/cff (</>), 



(5) 




FIG. 1: Phase-coupled flux qubits with a connecting loop in- 
terrupted by two dc-SQUID's. The arrows indicate the flow 
of the Cooper pairs and thus in reverse direction is the cur- 
rent. Gray squares denote Josephson junctions with Joseph- 
son coupling energy Eji for qubit loops and E'j for connecting 
loop. The qubit operating point is ft ~ 0.5 and fn ~ 0.5. 
(a) Two phase-coupled qubits are in the same current state, 
| ||), where two phases <f>tl and cj>m are nearly equal to each 
other resulting negligible coupling energy. Here | |) and | |) 
denote the diamagnetic and paramagnetic current state, re- 
spectively, (b) Left (right) qubit is in the diamagnetic (para- 
magnetic) current state, | ||), where the large phase differ- 
ence, 4>l\ — <t>m, induces large Josephson energy and current 
in the connecting loop. 



which describes dynamics of a particle with effective 
mass M in the effective potential U e s{4>) with <f> = 
(<j>Li, <t>L2,4>L3,4>Ri,4>R2, 4>R3, <f>[, <t>2, <t>a, 4>a)- The kinetic 
part of the Hamiltonian comes from the charging energy 
of the Josephson junctions such as 



^) = 5(S) £ £ C «& +C S4> .(6) 

P=L,R \i=l / 



where Cnn)i and C is the capacitance of the Josephson 
junctions of the left (right) qubit loop and the connecting 
loop, respectively. The number of excess Cooper pair 
charges on Josephson junction iV, = Qi/q c is conjugate 
to the phase difference fa such as [fa,Ni] = i, where 
Qi = Ci( < I ) o/27r)0i, q c = 2e and Ci the capacitance of 
the Josephson junctions. Here we introduce the canonical 
momentum Pj and the effective mass My 

Pi = iV.fi = -ifi— Mij = 77^ CkSij, (7) 
dfa V 27r / 

to obtain the kinetic part of the Hamiltonian. 

The effective potential of the coupled qubits is com- 
posed of the inductive energy of loops and the Josephson 
junction energy terms; 



U eS (cf>) = U iad (<f>) + C7 qubit (0) + U conn (4>), (8) 

(9) 



U ind {<t>) = \Ls{ll + ll) + \L'I' 2 



Uqnbit(4 > )=y^E Ji (l - COS(f> Li )+^2Eji(l - COS^fli) 
i=l i=l 

(10) 
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cw*)=yX(i-cos#)- (ii) 



z=l 



Here U- m d(<f>) is the inductive energy of loops with the 
current of the right qubit Ir, left qubit II and connect- 
ing loop /'. f/ qu bit(0) is the energy of the Josephson 
junctions in two qubit loop and U conn (4>) that of the con- 
necting loop with Josephson coupling energies Eji and 

E'j- 

In experiments the two Josephson junctions with phase 
differences 4>l(r)2 and 4>l(r)3 can be considered nomi- 
nally the same so that it is reasonable to set 

Ej2 = Ej a = Ej, 4>L(R)2 = 4>L(R)3- (12) 

Here we introduce a rotated coordinate 



4>p = (0L3 + <Ar3)/2, 
fan = (4>L3 ~ (l>R3)/2 



(13) 

(14) 



and then using the boundary conditions in Eq. Q wc 
get 



?L1 



± (/>m = -40 p(m) + 2?r(n L ± n R + f L ± f R ). (15) 
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Thus we can reexpress the sum of Josephson junction 
energies of both qubits as 

£/qubit(</>) - 2E J1 [1 - C0S(P7T - 2<j) p ) cos(Af7r - 2(j> m )} 

+4Ej(l — cos (p p cos0 m ), (16) 

where P = n L + n R + f L + fn + fmd,L + /md,fl and M = 
n L - n R + f L - f R + /ind,L - /ind,fl- Since experimentally 
qubit operations are performed at near the co-resonance 
point fi, = f R = 0.5 and the induced flux is so weak 
as /ind,L(R) ~ 0.002, P and M can be approximated as 
integers such that P = + hr + 1 and M = til — n R . If 
P is even, M is odd and vise versa, so we can get simple 
form for [/ qu bit(0), 



U, 



qubit \}Pm > (Pp 



) = 2Ej\ cos 2<j>p cos 2cj) r , 



—4:Ej cos 4>p cos 



2Eji + 4Ej. 



(17) 



(18) 
(19) 



Introducing another rotated coordinate 
<t>' P = ifti + 3 )A 

4>'m = (<t>'l - 03)/ 2 

and using the boundary conditions in Eq. Q to get 

± &)/2 = ^ (m) + n(f L T f R+PTq ), (20) 

the Josephson junction energy of the connecting loop 
U colm ((f>) = 53i=i — cos<^) can also be written as 

U coan {4>' m ,<f)' p )/2E'j = 2 - cos 0^ cos 0; (21) 
- co S [ft m + n(f L + f' R )} cos[^ + n(f L - f' R )], 

where we set p = and q = 0. 



III. COUPLED-QUBIT STATES IN EFFECTIVE 
POTENTIAL 

First of all we consider the case that the control fluxes, 
f' L and f' R , have opposite directions such that 



}l — Ir — f ■ 



(22) 



Note that the boundary conditions in Eq. ® already 
have opposite signs. In order to obtain the effective po- 
tential as a function of (f> m and <f> p , we reexpress $ as 
ftp = n (fcs + r + /4d) = _2 ^m + irM' using the bound- 
ary conditions in Eq. J5J and the expression in Eq. I)15[l. 
Here M' = M+r+/j ad can be written as M' = nL—n R +r 
neglecting small induced flux / nld in the connecting loop. 
Depending on whether M' is even or odd, the results will 
be quantitatively different, but qualitatively the same. 
Here and after, thus, we choose M' is even and specifi- 
cally til — 0, nji = and r = for simplicity and then 



i)' p becomes 




0.5 0.5 

<J)' /2n 



FIG. 2: (Color online) Effective potential of the coupled 
qubits in Eq. (25) for E'j = Q.lEj when f L =f' R = f' = (a) 
in (<f> m ,(j>p) plane and (b) in (4>' m ,(f)p) plane. Coupled qubit 
states at the local minima of potentials are denoted in pseudo- 
spin notation, which shows that these states are stable in both 
planes. Here we set Eji = Ej and f L =f R = 0.5. 



and the energy of Josephson junctions in connecting loop 
in Eq. (23J 

U CO nn(<fi' m , <p m ) 

= AE'j[l - cos7r/'cos(^ n + tt/') cos20 m ]. (24) 

Since the induced energy U- m d{4>) can be negligible, the 
total effective potential U c r(4>) in Eq. © is given by the 
sum of the energies in Eqs. I|17fl an d (|24|l . 

U e e(4>m,4>p, (f>'m) = ^qubit(^ )■ (25) 

The lowest energy level of U e g in (</> m , 4> p ) plane can be 
obtained by setting the remaining variable 4>' m in Eq. I|24|) 

as 

(i) 4>' m = 0, for - tt/4 < 4> m < tt/4, (26) 

(ii) 4>' m = ±7T, for tt/4 < \4> m \ <tt/2, (27) 

when /' = 0. We plot the effective potential U c s(4> m , <f> p ) 
in Fig. 12 (a) with four local minima. 

The value of local minima of case (i) can be obtained 
from U e s{4> m = 0,<fi p ) and we have found that two local 
minima have the same value, E ss (f), for equal pseudo- 
spin state with s S {j, f }. Similarly we get E s> - S (f) of 
case (ii) from U e s(4> m , 4> p — 0) for different pseudo-spin 
state. As a result, we obtain 

E ss ,(f) = iEj + iE'jil-cosnf) 

- 2Ejcosct) ss ,, (28) 

where 4> ss is the value of (j) p at local minima of the same 
spin states, \ss), and 4> s ,-s the value of 4> m of the different 
spin states, \s, —s) with 



E, 



(23) 



COS( 



2(^j! +2£^ costt/') 



(29) 
(30) 



4 



Thus the energy of the same spin states, E ss (f = Q), is 
lower than that of different spin states, E s _ a (/' = 0), as 
shown in Fig. Eta). 

Here we set Ej\ = Ej, E'j = O.lEj and f L = f R = 
0.5. For the same spin states we have two solutions, 
fa s /27r = ±1/6, corresponding to two local minima, 
E n (f = 0) and E n (f = 0). When fa s /2TT = 1/6, 
fa/2ir =1/6 and <p m — and thus 

<\>L(R)ll^ = <f>L{R)3/^ = <t>L{R)l/^ = V 6 ( 31 ) 

using 4>l(r)i + %4>l(r)z = 7r from the boundary condition 
in Eq. JlJ with num = 0. Since the loop currents of 
both qubits then 

I = -{2ir/$ )Ejsm4> L{m , (32) 

are diamagnetic as can be seen from Fig. [fla), this 
coupled-qubits state can be represented as | J.J.) as shown 
in Fig. Efa)- On the other hand, when 4> ss /2tt = —1/6, 

0L(fl)2/27r = <p L{m /2ir = -1/6, (f> L (R)i/2ir = 5/6. 

(33) 

Then the qubit current I = — (2-k/$o)Ej sin(— 7r/3) = 
— (2tt /$o)Ej sin(57r/3) corresponds to the paramagnetic 
current states, | TT)- We would like to note that, since 
the external fluxes /z, and f R threading left and right 
qubit loops are already in opposite directions, diamag- 
netic (paramagnetic) currents of both qubits in | J. J.) 
(| |T)) state are also in opposite directions. 

For different spin states, two solutions are also ob- 
tained for fa- s /2n w ±0.181. When fa ,- s /2tt « 0.181, 

fa 2 /2TT = (t> L3 /2wn 0.181, </> L i/2n « 0.138 (34) 
<f> R2 /2n = c/) R3 /2tt w -0.181, (f> R i/2w w 0.862 (35) 

for left and right qubit respectively, which corresponds 
to the state, | ||) in Fig. El ( a ). I n the same way 
4>s.- s /2tt « —0.181 corresponds to the state | Tl). Hence 
we can identify four stable states, | ||), | ff), | ||) and 
Tl), with energies E ss and E S: _ S at the local minima of 
U eS (4> m , fa) as shown in Fig. El (a). 

Even though above four states are stable states in 
(</>,„, </> p ) plane, it can be unstable in the other dimen- 
sions if they are saddle points. Thus we need represent 
the effective potential U c g in [fa m ,(f)') plane. From the 



\ss') 




(</> m /27r, 4> P /2n) 


(0^/2jr,#,/27r) 


1 II) 


(1/6,1/6) 


(0, 1/6) 


(0,0) 


1 TT) 


(-1/6,-1/6) 


(0,-1/6) 


(0,0) 


1 IT) 


(-0.181,0.181) 


(0.181,0) 


(±0.5, -0.362) 


Tl) 


(0.181,-0.181) 


(-0.181,0) 


(±0.5,0.362) 



TABLE I: The values of phase differences of coupled qubits 
states in several coordinates with E,j\ — E.j, E'j — O.lEj, 
h = f R = 0.5 and /' = 0. 



expression of U c s((f> m ,(t)p,(j)' m ) in Eq. (|2l)|l and the rela- 
tion <j>' p = -2<f) m in Eq. ESI, we can get U cS {(j)' m , <j/ p , fa) 
and, following similar procedure as in the (4> m , fa) plane, 
we obtain the effective potential as shown in Fig. ElO 3 ); 
where we can again see local minima. In Figs. El(b), for 
the states | ||) and | TT) of case (i), we can get the values 

fa p = and fa m = 0, (36) 

and, for case (ii), the values 

f fa p /2ir = -0.362, fa m /2ir = ±0.5, for | ||) , 
1 0;/2tt = 0.362, ^/27r = ±0.5, for | Tl) 

using Eqs. (|23|l . 126|l and H27|l and the values of 4> m in 
each case. As a result, we are able to identify the spin 
states at local minima of Figs. El(b) from Fig. El( a ) with 
above values and confirm the stability of the states in 
both planes. In Table U we summarize the values of the 
phase differences for four states, \ss'), of coupled qubits 
in several coordinates. Actually we obtained higher en- 
ergy states in Fig. El(a), but found that they are unstable 
in P lane - 

IV. TUNABLE COUPLING OF FLUX QUBITS 

The Hamiltonian of coupled qubits can be written as 

#cou P = h L a z L <g)I + h R I ®a R - Ja z L ® a z R 

+ t L a x L ® I + t R I ®<t r + E , (38) 

where h L = {E^+E^)/2-Eq and h R = {E ll +E 11 )/2- 
E with E = (Etf + E n + E Ll + E n )/4 and / is the 
2x2 identity matrix. First two terms are qubit terms, 
the third is coupling term and last two terms are tun- 
nelling terms which come from the quantum fluctuation 
described by the kinetic term of the Hamiltonian. Then 
the coupling constant J of the coupled qubits is given 



In Fig. |3 we plot the energies of coupled-qubits for 
various /' with E'j = O.lEj, Eji — Ej and /x = f R — 
0.5 in {4>m,<j>Lz) plane. When /' = in Fig. Eta), the 
energies E ss of the same spin states, | ||) and | TT)> arc 
lower than E s _ s , of the different spin states, | ||) and 

Tl)- The positions of four local minima are shown in 
Table 0] As increases /', the energy difference AE = 
i? s ._ s — E ss becomes smaller (upper panel in (b)) and 
finally AE = at /' = 0.5 in (c). Since E u = E n = E ss 
and Eii = ^Ti = ^s.-s, the coupling strength can be 
written as 

2J(f) = AE(f') = E Si _ s (f) - E ss (f). (40) 

Therefore the coupling strength between two flux qubits 
changes as varying the control fluxes /' threading the 
dc-SQUID loop in the connecting loop. 
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-0.2 0.2 

(b) 



FIG. 3: (Color online) Effective potential of the coupled 
qubits in Eq. © for E'j = 0.1K/, Eji = Ej and ft = 
f R — 0.5. Coupled qubit states at the local minima of poten- 
tials are denoted in pseudo-spin notation, (a) Effective po- 
tential as a function of <f) R 3 and <j>t3 when f' L = f' R = /' = 
for the phase-coupled qubits in Fig. Here the energies of 
different current states are equal to each other, E^ = Efi, 
as well as En — E^ for the same current states. The en- 
ergy of different current states E s ^ s = Ey\ = E-\y is higher 
than that of the same current states E sa = E^ — E^f. (b) 
(top) Two control fluxes in Fig. Qare in opposite directions, 
ft = /r = /', an d /' is increased to /' = 0.25. The energy 
difference AE — -E s ,_ s — E as becomes smaller than when 
/' = in (a), (bottom) For the case when two control fluxes 
are in the same direction such as f' L = — f' R = /' = 0.25, the 
energies of different current states are not equal to each other 
any more; E^ > E^. (c) The coupling becomes switched off 
when f' L = = /' = 0.5. Thus the energies of four states 
have the same value, En = En = E^ = En. 




-0.50 -0.25 0.00 0.25 0.50 



f 



-0.50 -0.25 0.00 0.25 0.50 



f' 



FIG. 4: (Color online) Energies of coupled qubit states for 
E'j = O.lEj, Eji — Ej and ft = f R — 0.5 as a function 
of /'. (a) E ss ,(f) in Eq. 128H when two control fluxes in 
Fig. 0are in opposite directions, f' L = f' R = /'. The coupled 
qubits can be described by the Hamiltonian in Eq. 1381 and 
the coupling strength 2J(/') = E a - a {f) - E ss (f) in Eq. 
I|41|l is also shown. As increases /', the coupling strength 
decreases monotonously, vanishing finally at /' = 0.5. (b) 
When two control fluxes in Fig. 0are in the same direction, 
ft = -fn = /', the energy of two states, | J.T) and | ||), 
becomes different, as if additional fluxes, 5ft and 5f R , were 
applied into the left and right qubit loop, respectively. Hence 
the coupling between two qubits cannot be represented solely 
by change of the coupling constant of the Hamiltonian in Eq. 
l|38jl . (c) For the case when there is only one dc-SQUID loop 
in the connecting loop instead of two dc-SQUID's in Fig. 
the energies of the different current states, | jf) and | |J.), are 
also different from each other. 



From Eq. Q28JI the coupling constant J can be rep- 
resented as a function of /' by J(f') — Ej(cos(f> ss — 
cos 4> s which gives 



J(f') = 



Ej 



E'j COS IT f 



Ejx E J1 +2E'cosirf r 



(41) 



In Fig. 01a) we plot the energies E ss i{f) and J(f') as 
a function of /', where 2J(f) = E Si _ s {f) - E ss (f). 
When /' = 0, J is of the order of E'j so that we can 
obtain a sufficiently strong coupling. By adjusting /' the 
coupling strength can be tuned from strong coupling to 
zero at /' = 0.5. 

The coupling strength J(f') in Eq. I|41|) depends on 
Ej/Eji as well as E'j. When E'j/Ejx is small, J(f') is 
proportional to E'j and (Ej / Eji) 2 . Recently the phase- 
coupling scheme has been experimentally implemented^, 
where four-Josephson junctions qubits are employed in- 
stead of usual three- Josephson junctions qubits. In that 
experiment the Josephson junction energy Ej\ of fourth 
junction is large so that the value of Ej/Eji is about 
Ej/Eji w 1/3. As a result, the experiment exhibits 
rather small coupling strength. 

The current of connecting loop can be written as 



-($o/27r)/' = ^£i=i sin $ - ^.'/Eta 8111 ^' which 
gives the relations, (^ — ^3) + {(t>' 2 — 4>'i) — 47rfc and then 



<P' m = (k- />, 



(42) 



with integer k using the boundary conditions in Eq. 
©. Then, using Eq. J20I) and the effective flux f' cS 
in Eq. (23, the current -($ /27r)/' = Q.5E'j Y% =1 sin$ 
is given by 



f = (19 (-^^cosTr/'sin^. 



(43) 



This current-phase relation can be consid- 
ered as the Josephson junction type relation, 
V = -(27r/$o)(-l) fc -E,/Sin^, with the effective 
Josephson coupling energy, E'j, of two dc-SQUID's in 
the connecting loop 



E'j = 2E'j cos nf 



(44) 



and the phase difference if = ftf^- The coupling con- 
stant in Eq. I|41l) also can be represented by the effective 
Josephson coupling energy, E'j. Thus the large phase 
difference, 7r/^ ff , and the Josephson coupling energy, E'j, 
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induce the current in the connecting loop and the cou- 
pling energy of the phase-coupled qubits. 

For the same spin states, | J, J,} and | TT)- the current 
of connecting loop I' becomes zero, since 4>n — 4>ri 
and thus f' eS — 0. For a different spin states | If) 
with /' = 0, f eS = (4> tl ~ <t>Ri)/2ir » -0.724 and 
we have k = 1 from <fi' m = tt and the relation in Eq. 
(|42[1 . Then weak current I' in the connecting loop flows 
satisfying current conservation condition between left 
qubit and connecting loop such that Ej sin0.181(27r) = 
Ejt sin0.138(27r) + 2^sin0.724?r for /' = 0. When /' 
approaches 0.5, the effective Josephson coupling energy 
Ej = E'j cos 7r/' and thus the current I' in connecting 
loop become zero, which means that the coupling be- 
tween two qubits is switched off. 

Now we want to explain the case that two control fluxes 
are in the same directions and the case that there is only 
single dc-SQUID in connecting loop. If two control fluxes 
are in the same direction such as 

S'l = S'r = (45) 

the Josephson junction energy of the connecting loop be- 
comes 

U conn {^' m ,(j) m )=^E'j[l - cos tt/' cos 4> m cos(20 m -nf')]. 

(46) 

Similar procedure as in the case of opposite directions of 
control fluxes shows that the same spin states, | J, J.) and 
| It), have equal energy such as 

E U = £ TT = Ej ^4 - + 4E'j sin 2 tt/'. (47) 

for cos (j) p — Ejj1Ej\. 

For different spin states, | if) and | fj,), the energies 
En ano - -^it are obtained at two local minima 

U conn (4> m ,(t>p = 0) = -iE'j cos 2 tt/' cos 2(j) m 

-2E'j sin 2tt/' sin 2<j> m + AE'j, (48) 

which can be derived from Eq. (|46[) . Since the states, 
I It) an d I Tl), have different sign for <j)s,-s, the second 
term produces the energy difference 

AE = E n - E ix = 4S' 7 sin27r/'| sin2^ s _ a \, (49) 

where Si _ s is again one of the values of <j) m for the dif- 
ferent spin states. 

Figure 01b) (lower panel) for /' = 0.25 shows that, 
when two control fluxes are in the same direction, the 
energies Ef± and E^ are different while En = E^. The 
energy levels of E ss > are plotted in Fig. Efb). In this 
case the effective fluxes hi and h R applied to left and 
right qubits in the Hamiltonian of Eq. I|38[) become dif- 
ferent each other, hi ^ h Rl as /' increases from zero. 
For the different current state in Fig. ^b) , if the control 
fluxes f' L and f' R threading the dc-SQUID loops are in the 



same direction, the increased current I' in the connect- 
ing loop will flow through the left and right qubit loops. 
Thus the qubit states are influenced by additional effec- 
tive fluxes, which will makes the two-qubit operations 
difficult. However, if two control fluxes f' L and f' R are in 
opposite directions, the energies of different spin states 
remains equal to each other, E^i — E^, as shown in Fig. 
St a). This means that the additional currents coming 
from two dc-SQUID's are cancelled each other and total 
additional current induced by the control fluxes f' L and 
f' R is vanishing in the connecting loop. As a result, the 
net effect is just renormalizing the coupling constant J 
of the coupled qubit system. 

We also calculated energies of coupled qubit states 
with single dc-SQUID loop whose boundary conditions 
become 

0i = 27r(r + /4 d ) + (^n-^m) (50) 
-<P' 1 + <t>' 2 = 2Tr(f'+p), (51) 

instead of those in Eqs. and Then we get the 
Josephson junction energies of the dc-SQUID, 

U conn (4> m ) = -2E'j cos(4<£ m - tt/') + 2E'j, (52) 

which gives results similar to those of two dc-SQUID's 
with fluxes in the same direction such that 

E u = E n =Ej(a- + 2E[ } sin 2 tt/' (53) 

for cos (j)p — Ej/2Eji and 

AE = E n - E u = 2^sin27r/'|sin40 Si _ s | (54) 

as shown in Fig. 0fc). Hence the behaviors of one dc- 
SQUID in the connecting loop are qualitatively the same 
as those of two dc-SQUID's with fluxes in the same direc- 
tion. Therefore we need two control fluxes threading dc- 
SQUID's in opposite directions to cancel the additional 
currents in the connecting loop for obtaining the control- 
lable coupling. 

In order to obtain the controllable coupling both the 
qubit operating flux, fi, and control flux, f' L , of the left 
qubit become in opposite direction to those of the right 
qubit, f R and f' R as shown in Fig. ^ In real experiments 
it will be very hard to apply magnetic fluxes of different 
directions simultaneously. We have previously suggested 
a scalable design for phase-coupled flux qubits^, where 
an arbitrary pair of qubits are coupled in a twisted way. 
Thus just applying all magnetic fluxes in the same direc- 
tion makes automatically the effect of fluxes in opposite 
directions, removing the experimental difficulty. 

The recent experiment on the phase-coupled flux 
qubits without dc-SQUID loopA& has shown that the cou- 
pled qubit states are in quantum mechanically super- 
posed regime. The dc-SQUID loops in the connecting 
loop of the present tunable coupling scheme may cause 
a decoherence effect on the coupled qubit states. A re- 
cent study argued that the dc-SQUID based oscillator 
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should be the main source of the decoherence of the flux 
qubits^. For the scalable design in Ref. 0, however, the 
decoherence from two dc-SQUID's can be reduced. Since 
two dc-SQUID's are connected in a twisted way, the fluc- 
tuations from tank circuit or flux lines can be cancelled 
each other. 

In realistic implementation of qubit operations, oper- 
ating external fluxes are slightly different from the co- 
resonance point, fh = fn = 0.5, and moreover we cannot 
any more neglect small kinetic inductance and induced 
fluxes. Hence we confirmed the results in this study by 
numerical calculation using the exact boundary condi- 
tions similar to those in Eqs. Q — ©, current-phase re- 
lation Ii = — (2n/$Q)Eji sm<j)i and current conservation 
conditions^. 



V. SUMMARY 

Controllable coupling between two phase-coupled flux 
qubits can be achieved by using two dc-SQUID's in the 



connecting loop with threading fluxes in opposite di- 
rections. We analytically show at co-resonance point 
(/l = Ir = 0.5) that the coupling strength of the 
phase-coupled flux qubits can be adjusted by varying the 
threading fluxes /' from to 0.5; it can be as strong 
as O(E'j) and zero in switching-off limit. When cither 
two control fluxes are in the same directions or there is 
only one dc-SQUID in the connecting loop, the coupled 
qubits cannot be described by the coupling Hamiltonian. 
In slightly different parameter regimes of experimental 
implementations numerical calculations can be done to 
obtain exact results. 
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